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Abstract: An observation of seasonal variations in underground muon rate, i? M , has been performed at 
Soudan, MN, by the MINOS Far Detector. The four percent fluctuation seen over three years was highly 
correlated to the temperature variations of the upper atmosphere. The coefficient relating variations in 
temperature to variations in muon rate was found to be: cut = (T / R l± )(dR ll /dT) = 0.87 ± 0.03, which 
is near the expectation of 0.91. 



Introduction & Motivation 

When cosmic rays interact with molecules in the 
Troposphere, mesons are produced which either 
interact again and produce low energy cascades 
or decay into muons. While the temperature of 
the Stratosphere varies considerably within the day 
in areas far from the equator, the temperature of 
the Troposphere remains nearly constant, slowly 
changing over longer timescales such as seasons. 
Increases in temperature of the Troposphere cause 
increases in volume and atmospheric scale height, 
thus the height of the primary cosmic ray interac- 
tion. The higher in the atmosphere mesons are pro- 
duced, the more time they have to decay to muons, 
thus the rate of muons underground will increase 
as temperature increases [1-3]. Though this effect 
has been measured by underground experiments, 
there has been little agreement with the expecta- 
tion. 

MINOS is a long baseline neutrino oscillation ex- 
periment, with a beam and Near Detector at 
Fermi National Accelerator Laboratory in Batavia, 
IL. The Far Detector is a 5.4 kt magnetized scintil- 
lator and steel tracking calorimeter located 720 m 
underground (2100 mwe) at the Soudan Under- 
ground Mine State Park in Northern Minnesota. Its 
depth, large acceptance and flat overburden make 
it possible to observe cosmic -ray induced muons 
of minimum surface energy 0.7 TeV without pref- 
erence to direction, and thus detect the small sea- 



sonal fluctuations in arrival rate. The seasonal ef- 
fect is enhanced as muon energy increases, and 
the large size of the Far Detector allows a signif- 
icant accumulation of statistics with which to per- 
form this analysis. Additionally, the Far Detector 
has a magnetic field, which allows the separation 
of particles by charge, so MINOS will be the first 
experiment to measure seasonal variations for /i + 
separate from The consistency and availabil- 
ity of radiosonde temperature measurements from 
the NOAA IGRA (Integrated Global Radiosonde 
Archive [4]) over the duration of the data set en- 
sures a high statistics temperature sample as well, 
increasing the probability of a positive correlation. 
The data used in this analysis were collected over 
three years, from 1 August 2003 - 1 August 2006 
for three complete cycles, numbering 24 million 
muons. The relationship between the temperature 
and intensity can be expressed as [1]: 

AI U f°° , ,AT(X) 

where AI^ are the fluctuations about 7°. The short 
lived mesons produced in the upper atmosphere in- 
teract or decay as they descend toward the earth. 
The meson decay channels result in muons with 
nearly the same energy as the parent meson, while 
interactions produce lower energy cascades that are 
filtered by the rock overburden above the Far De- 
tector. These outcomes are energy dependent, sep- 
arated by the "critical energy" [2]. The "Effec- 



Observation of Seasonal Variations. 



tive Temperature" (T e //) approximates the upper 
atmosphere as an isothermal body, weighting the 
temperature of the pressure levels to have a uni- 
form amount of matter. In the ir scaling limit, T e f t 
is [2]: 



T, 



e// 



(2) 



where X is the scale height of the atmosphere, 
An — 120gm/cm 2 and A w = 160gm/cm 2 
are the nucleon and pion atmospheric attenua- 
tion lengths, respectively. For a detector count- 
ing discrete particles, the intensity is written I ^ — 
Ri/eA e f ffl, where Ri — Ni/ti, the number of 
muons observed over time ti,A e ff is the effective 
area, e is the efficiency, and ft is the solid angle 
observed. Every term but the rate is constant over 
time, so: ~ jr^x. With these definitions and 
eq. 1, we can write the experimental determination 

of (XT'- 



The Data 
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The data for this analysis were accumulated over a 
three year span, beginning on 1 August, 2003, at a 
time when the detector was fully operational. Be- 
ginning with 40.3 million cosmic ray tracks, a se- 
ries of cuts were performed [5]. Pre-analysis cuts 
include: failure of demultiplexing figure of merit, 
multiple muon (multiple muons aren't included in 
the Monte Carlo), "bad run" and bad magnet coil 
status. Analysis cuts include: track length less than 
2 m number of planes less than 20, Xreco > 1 and 
either track vertex or end point outside of the fidu- 
cial volume of the detector. A total of 24.7 million 
events survived these cuts for the combined sam- 
ple. T e f f was found using weather data from Inter- 
national Falls, MN weather station. Balloon flights 
were usually done twice a day, with the maximum 
height reached at noon and midnight, and sampled 
temperatures from at least six different heights. 
Days in which there were not exactly two temper- 
ature readings or that both measurements did not 
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Figure 1 : Time between consecutive [i arrivals in 
log y. A Poisson fit gives x 2 /ndof = 90.87/105; 
(R,j,) (from slope) = 0.2872 ± 0.0001 
reach a column depth of at least 60gm/cm 2 were 
excluded from the data set. 



Analysis 

Upon examination of the data, it was found that on 
four days there were fluctuations that deviated in 
an erratic manner. The great stability of the detec- 
tor over the 1096 days of data and the fact that they 
were documented by the Control Room Logbook 
made these days stand out and diagnose as hard- 
ware issues. To find the rate for each day, the num- 
ber of muons counted was divided by that day's 
livetime. T e f t was calculated for two times each 
day using the IGRA temperature data and 2, and 

the error was found by er 2 = (t^^ — (T e //) 2 
added in quadrature with 0.05°. The fit results 
from Fig. 1 was used to find (Ru) over three years, 
0.287 Hz. Histograms of the deviations from the 
mean for both R^ and T e f t are shown in Fig. 2, 
binned by day. The expected periodic fluctuation 
in T e ft, with maxima in July, minima in January, 
is very clearly shown, as is a very similar (nearly 
indistinguishable) fluctuation in R^. An indepen- 
dent analysis used a smoothed time series, and 
their results were highly consistent with what is 
shown here. To quantify the correlation between 
rate and temperature, a plot of R^(T e ff) was pro- 
duced (Fig. 3) and a linear regression was fit using 
ROOT'S MINUIT fitting package. This package 
accounts for error bars on both the x and y axis 
using a numerical minimization method. This fit 
gives a-r = 0.87 ± 0.03 from the slope. In order 
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Figure 4: The theoretical cxt{X) (solid curve) for slant depths up to 4000 mwe. The MINOS point is from 
this analysis, Barrett 1,2 [1], AMANDA [3]; all other points from [2] 
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Figure 2: Ai? p , (top), and AT e // (bottom) from 
8/03-8/06, binned by day. 

to compare our experimental cut to the theoretical 
expectation, a simple numerical program was writ- 
ten to find the expected value given by [ 1 ] : 



(«t) w = (1 



7 



(4) 



( 7 + l) l.\E th cos0 

Note that this expression is only valid for pions. 
Future work will involve this kaon contribution, 
which should lower the expected olt since kaons 
are short lived and always decay. A muon energy 
and cos were chosen out of the differential muon 
intensity [6], 
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Figure 3: A plot of Ai? M / vs. 
AT e ff/ (T e ff) for single muons. The fit 
X 2 /ndof = 1420/953, and correlation coefficient 
R = 0.79. 

where 7 = 1.7 is the muons spectral index [5]. a 
random azimuthal angle, <f> was chosen and com- 
bined with cos 8 and Soudan rock overburden map 
[7] to find the slant depth. The threshold surface 
energy required for a muon to survive this column 
depth is found from Eth(&, 4>) = a(e hx — \), where 
a = 0.45 TeV and b = 0.44 [kmwej^ 1 for Soudan 
rock [5], column depth X = X(9, </>), and if the 
chosen E^ > E t h, it was used in the calculation 
of the theoretical (ar)^- This was repeated for 
10,000 successful E^ to find (a^ = 0.91 for 
MINOS, which is very near to the experimental 
value, 0.87 ± 0.03. To compare the MINOS re- 
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suit with other underground experiments, this pro- 
cess was repeated for standard rock (a = 0.50 TeV 
and b = 0.4 [kmwe] -1 ), fiat overburden, and X 
II cos 8, where H is the detector depth in mwe, 
using 10,000 successful muons at depths from to 
4,000 mwe. The result of this calculation, along 
with data from other experiments, can be seen in 
Fig. 4. The MINOS result matches the expectation 
and has tighter error bars than both recent results, 
AMANDA (±0.05 [3]) and MACRO (±0.13 [2]). 

The curvature of the track is used to determine 
the momentum and charge of the particle, so a 
charge sign confidence cut was required. This cut 
was charge over momentum divided by the error 
in the determination of charge over momentum 
( MlL > 2.2), determined from previous investiga- 
tions of the muon charge ratio. That left 8.8 million 




Figure 5: AR^/ for fi + (open triangles, top) 
and /i~ (open circles, bottom), binned by day. 

events; 5.1 million positive, 3.7 million negative, 
which is consistent with the published MINOS 
charge ratio. Fig. 5(t) shows Ai? M + (open trian- 
gles) and Fig. 5(b) shows Ai? M + (open circles) over 
the same time period, binned by day. The sample 
of muons is smaller than for the fi to t sample, thus 
the error bars on are larger, but the trade off 
is that the error bars on the temperature are much 
smaller since the small fluctuations over several 
days are not washed out. Performing the same fit 
of Rfi(T e ff) as for the /itot sample on /i + sepa- 
rate from /i~ resulted in a slope of 0.845 ± 0.036 
and 0.843 ± 0.042 respectively. These correspond 
highly to each other, and are within one sigma of 
ax found from the pL to t sample. 



Conclusions 

A three year sample of 42 million cosmic ray in- 
duced muons has been collected by the MINOS 
Far Detector and daily rate fluctuations have been 
compared to daily fluctuations in atmospheric tem- 
perature, and these distributions were shown to be 
highly correlated, with a correlation coefficient of 
0.79. The constant of proportionality relating the 
two distributions, ax, was found to be 0.87 ±0.03, 
which, within the error band, is in good agreement 
with the theoretical expectation in the pion-only 
approximation of (ar)^ — 0.91. This suggests 
that the majority of muons seen in the Far Detector 
were generated by pion parents. 
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